BLACK HOLES, GALAXY FORMATION, AND THE Mbu - a 

RELATION 



Andrew King 



1,2 



ABSTRACT 

Recent X-ray observations of intense high-speed outflows in quasars suggest 
that supercritical accretion on to the central black hole may have an important 
effect on a host galaxy. I revisit some ideas of Silk and Rees, and assume such 
flows occur in the final stages of building up the black hole mass. It is now possible 
to model explicitly the interaction between the outflow and the host galaxy. 
This is found to resemble a momentum-driven stellar wind bubble, implying a 
relation Mbh = (/g't/27rG'^)cr^ ~ 1.5 x lO^cigoo Mq between black hole mass and 
bulge velocity dispersion {fg = gas fraction of total matter density, k = electron 
scattering opacity), without free parameters. This is remarkably close to the 
observed relation in both slope and normalization. 

This result suggests that the central black holes in galaxies gain most of their 
mass in phases of super-Eddington accretion, which are presumably obscured 
or at high redshift. Observed super-Eddington quasars are apparently late in 
growing their black hole masses. 

Subject headings: accretion - quasars: general - galaxies: formation, nuclei - 
black hole physics 



Introduction 



It is now widely accepted that the centre of every galaxy contains a supermassive black 
hole. The close observational correlation (Ferrarese & Merritt, 2000; Gebhardt et al., 2000; 
Tremaine et al., 2002) between the mass M of this hole and the velocity dispersion a of the 
host bulge strongly suggests a connection between the formation of the black hole and of the 
galaxy itself. 
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Recent XMM-Newton observations of bright quasars (Pounds et al., 2003a, b; Reeves et 
al., 2003) may offer a clue to this connection. These observations give strong evidence for 
intense outflows from the nucleus, with mass rates Mout ~ IM© yr~^ and velocity v ~ 0.1c, 
in the form of blueshifted X-ray absorption lines. Simple theory shows that the outflows are 
probably optically thick to electron scattering, with a photosphere of ~ 100 Schwarzschild 
radii, and driven by continuum radiation pressure. In all cases the outflow velocity is close 
to the escape velocity from the scattering photosphere. As a result the outflow momentum 
flux is comparable to that in the Eddington-limited radiation field, i.e. 

MoutV ^ (1) 

where Afout is the mass outflow rate and L^dd the Eddington luminosity, while the mechanical 
energy flux is 

1 TiV 2 -^Edd /r,\ 

It appears that such outflows are a characteristic of super-Eddington accretion (King & 
Pounds, 2003). We know that most of the mass of the nuclear black holes is assembled by 
luminous accretion (Soltan 1982; Yu & Tremaine, 2002). It seems likely that the rate at 
which mass tries to flow in towards the central black hole in a galaxy is set by conditions far 
from the hole, for example by interactions or mergers with other galaxies. It is quite possible 
therefore that super-Eddington conditions prevail for most of the time that the central black 
hole mass is being built up. 

This clearly has important implications for the host galaxy. Unlike luminous energy, a large 
fraction of a mechanical energy flux like (2) is likely to be absorbed within the galaxy, and 
must have a major effect. To reach its present mass the black hole in PG1211+143 could 
have accreted at a rate comparable to its current one for ~ 5 x 10^ yr. During that time, an 
outflow like the observed one could have deposited almost 10^" erg in the host galaxy. This 
exceeds the binding energy ~ 10^^ erg of a bulge with 10^^ Mq and a ~ 300 km s~^. 

Accordingly it is appropriate to revisit some ideas presented by Silk and Rees (1998., hence- 
forth SR98) and also considered by Haehnelt et al. (1998), Blandford (1999) and Fabian 
(1999). These authors envisage a situation in which the initial black holes formed with 
masses ~ lO^M© before most of the stars. Accretion on to these black holes is assumed to 
produce outflow, which interacts with the surrounding gas. Without a detailed treatment 
of the outflow from a supcrcritically accreting black hole, SR98 used dimensional arguments 
to suggest a relation between M and a. However this still has a free parameter. Given the 
simple relation (1) one can now remove this freedom. The situation turns out to resemble a 
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momentum-driven stellar wind bubble. Modelling this gives an Mbh — c relation devoid of 
free parameters, and remarkably close to the observed relation. 



2. Black Hole Wind Bubbles 

I follow SR98 in modeUing a protogalaxy as an isothermal sphere of dark matter. If the gas 
fraction is fg — Jlbaryon /^matter — 0.16 (Spcrgcl ct al., 2003) its density is 

where a is assumed constant. The gas mass inside radius R is 

M{R) = 47r rpr^dr = ^^^-^ (4) 







G 



I assume that mass flows towards the central black hole at some supercritical rate Mace- The 
results of King & Pounds (2003) suggest that this will produce a quasi-spherical outflow 
with momentum flux given by (1). Note that this momentum rate is independent of the 
outflow rate Mout = ^acc — A^Edd since the outflow velocity v adjusts as to maintain 
the relation (1) (King & Pounds, 2003). 

The wind from the central black hole will sweep up the surrounding gas into a shell. As is 
well known from the theory of stellar wind bubbles (e.g. Lamers & Casinelli 1999) this shell 
is bounded by an inner shock where the wind velocity is thermalized, and an outer shock 
where the surroimding gas is heated and compressed by the wind. These two regions are 
separated by a contact discontinuity. The shell velocity depends on whether the shocked 
wind gas is able to cool ('momentum-driven' flow) or not ('energy-driven' flow). In the 
absence of a detailed treatment of a quasar wind, SR98 appear to have assumed the second 
case. In fact for the supercritical outflows envisaged here the flrst case is more likely, as the 
argument below shows. 



3. Cooling the Wind Shock 

The Compton cooling time of an electron of energy E is 



where me is the electron mass and 
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is the radiation density at distance R from the black hole, and b < 1 allows for some 
coUimation of the outflow. The electron energy E in the postshock wind gas is ~ dnipV^ 
where v is the wind velocity and nip the proton mass. Combining this with the usual 
definition 

AttGMbkc 

^Edd — [ I ) 

K 

of the Eddington luminosity for black hole mass Mbh shows that 

'-l^(:=;)'(D--=«^.©W'y. (s) 

where i?kpc is R measured in kpc and Mg = Mbh/IO^Mq. Clearly this is extremely short 
for small R, so the flow is efficiently cooled and thus momentum driven at least initially. I 
note that Ciotti & Ostriker (1997, 2001) emphasize the importance of Compton heating and 
cooling on quasar inflows and outflows. 

The momentum-driven assumption breaks down once tc becomes of order the flow time 
^flow — R/vs where Vs is the shell velocity. We can use the momentum-driven shell velocity 
derived in eqn (14) below to estimate 

tflow = 8 X lO^i^kpcaaooMg-'/' yr (9) 

where cr2oo = o'/(200 km s~^). The assumption of efficient cooling is valid out to a radius Rc 
given by setting tc — ^aow — 1- We find a total swept-up mass 

M(i?e) = 1.9 X 10^V^ooA4''^(^) M0 (10) 

at this point. Once the shell reaches radii larger than Rc the shocked wind is no longer 
efficiently cooled, and its thermal pressure accelerates the shell of swept-up gas to a higher 
velocity > Vm (energy-driven flow) after a sound-crossing time ~ Rc/v- 



4. The Mbh — cr relation 

I now estimate the speed Vm of the momentum-driven shell by the standard wind bubble 
argument. At sufficiently large radii R the swept-up shell mass M{R) is much larger than 
the wind mass, and the shell expands under the impinging wind ram pressure pv'^ (this 
characterizes momentum-driven flows; in an energy-driven flow the thermal pressure of the 
shocked wind gas is dominant, while in a supernova blast wave the momentum injection is 
instantaneous rather than continuous). The shell's equation of motion is thus 



d 

di 



M{R)R ^AnR^/w^ ^ M^^tv^^^ (11) 
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where we have used first the mass conservation equation for the quasar wind, and then (1) 
to simphfy the rhs. Integrating this equation for R with the final form of the rhs gives 

M{K)R=^t (12) 

where I have neglected the integration constant as M{R) is dominated by swept-up mass at 
large t. Using (4) for M{R) and integrating once more gives 

where again we may neglect the integration constant for large t. We see that in the snowplow 
phase the shell moves with constant velocity Vm = R/t, with 

We note that this velocity is larger for higher I/Edd, i-e. higher black hole mass. This solution 
holds if the shell is inside the cooling radius Re, outside this radius the shell speed eventually 
increases to the energy-driven value Vg, which also grows with Mbh- 

I now consider the growth of the black hole mass by accretion. Initially the mass is small, 
inflow is definitely supercritical, and even the energy-driven shell velocity would be smaller 
than the escape velocity o". No mass is driven away, and accretion at a rate MEdd can 
occur efficiently. However as the black hole grows, we eventually reach a situation in which 
Vg > a > Vm- Further growth is now only possible until the shell reaches i?c, and then only 
until the point where Vm = o. Thus given an adequate mass supply, e.g. through mergers, 
the final black hole mass is given by setting = c in (14). Using (7) we find the relation 

Mbh = ^^<7^ ^ 1.5 x lOV^^oo Mq. (15) 

This is remarkably close to the observed relation (Tremaine et al., 2000). 

Presumably most of the swept-up mass ends up as bulge stars, and we may tentatively 
identify M{R^ as an upper limit the bulge mass of the galaxy. Using (15) to eliminate 
c^2oo we get 

Mbh > 7 X lO-^Mg"'/' hM^. (16) 

If cjv (determined by the ratio Mout/MEdd) attains similar values at this point in most 
systems and the swept-up mass is close to M(i?c) one gets a relation between black hole 
and bulge mass of the form M(, cx Mg^^. The relation is written instead in the form (16) 
to allow easy comparison with the correlation found by Magorrian et al. (1998). Evidently 
this is not as clear-cut a relation as that between Mbh and cr, and indeed the scatter in the 
observed relation is considerably larger. 
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5. Discussion 

The Mbh — relation given here has no free parameter. If the outflow velocity v had been 
larger by an optical depth factor r > 1 (i.e. most of the acceleration occurs below the 
photosphere) a factor 1/r would appear on the rhs. However this would require outflow 
velocities t{GM / R^y^Y/"^ larger by the same factor than those observed in supercritically 
accreting quasars. 

The lack of freedom in (15) comes about because the physical situation envisaged by SR98 
and also studied by Haehnelt et al., 1998) and Blandford (1999) can now be made more 
precise: the response of observed black hole systems to super-Eddington accretion appears 
to be an optically thick outflow driven by continuum radiation pressure. Fabian (1999) 
considers su6-Eddington accretion, but emphasizes the importance of the momentum of the 
outflow as opposed to its energy: it is this which leads to the cr^ dependence rather than . 
Specifically, Fabian (1999) assumes a wind of speed with mechanical luminosity a fixed 
fraction a of LEdd- This produces a relation of the form (15) but with an extra factor tv./ac 
on the rhs; it therefore reduces to (15) if one assumes a ~ f^/c. Parameters also appear in 
other derivations using different physics, such as the ambient conditions in the host galaxy 
(Adams, Graff & Richstone, 2001) or accretion of coUisional dark matter (Ostriker, 2000). 

The picture presented here invokes a largely spherical geometry for the ambient gas, except 
that the accreting matter must possess a small amount of angular momentum to define an 
accretion disc plane and thus a small solid angle where inflow rather than outflow occurs. 
It is therefore appropriate to the growth of a spheroid-black hole system. However once 
most of the gas lies in the plane of the galaxy the momentum-driven outflow considered here 
would not halt inflow. Evidently this means that accretion from this point on adds little 
mass to the hole. 

If the derivation of the Mbh — cr relation given here is some approximation to reality, it 
implies that the central black holes in galaxies gain most of their mass in phases of super- 
Eddington inflow. As relatively few AGN are observed in such phases, these must either be 
obscured (cf Fabian, 1999) or at high redshift. It appears then that those quasars which 
are apparently now accreting at such rates (Pounds et al., 2003a,b; Reeves et al., 2003) are 
laggards in gaining mass. This idea agrees with the general picture that these objects - all 
narrow-line quasars - are super-Eddington because they have low black-hole masses, rather 
than unusually high mass inflow rates. 

I thank Andy Fabian, Jim Pringle and James Binney for illuminating discussions, and the 
referee for a very helpful report. This work was carried out at the Center for Astrophysics, 



- 7- 



and I thank members of its staff, particularly Pepi Fabbiano and Martin Elvis, for stimulating 
discussions and warm hospitality. I gratefully acknowledge a Royal Society Wolfson Research 
Merit Award. 

REFERENCES 
Adams, F.C., Graff, D.S., Richstone, D.O., 2001, ApJ 551, L31 

Blandford, R.D., 1999, in Galaxy Dynamics, ASP Conference Series 182, eds. D.R Merritt, 
M. Valluri, J. A. Sellwood (San Francisco: ASP), 87 

Ciotti, L., Ostriker, J.P., 1997, ApJ 487, L105 

Ciotti, L., Ostriker, J.P., 2001, ApJ 551, 131 

Fabian, A.C., 1999, MNRAS 308, 39 

Ferrarese, L., Merritt, D., 2000, ApJ 539, L9 

Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S.M., Filippenko, A..V., Green, 
R., Grillmair, C., Ho, L.C., Kormendy, J., Lauer, T.R., Magorrian, J., Pinkney, J., 
Richstone, D. and Tremaine, S., 2000, ApJ 539, L13 

Haehnelt, M., Natarajan, P., Rees, M.J., 1998, MNRAS 300, 817 

King, A.R., Pounds, K.A., 2003, MNRAS, in press (astro-ph/0305541) 

Lamers, H.J.G.L.M., CasineUi, J. P., 1999, Introduction to Stellar Winds, Cambridge Uni- 
versity Press, Cambridge, Ch. 12 

Magorrian et al., 1998, AJ 115, 2285 

Ostriker, J.P., 2000, Phys. Rev. Lett. 84, 5258 

Pounds, K.A., Reeves, J.N., King, A.R., Page, K.L., O'Brien, P.T., Turner, M.J.L., 2003, 
MNRAS in press (astro-ph/0303603) 

Pounds, K.A., King, A.R., Page, K.L., O'Brien, P.T., 2003, MNRAS, in press, astr- 
ph/0305571 

Reeves, J.N., O'Brien, P.T., Ward, M.J., 2003, ApJ, 593, 65 
Soltan, A., 1982, MNRAS, 200, 115 



-8- 



Spergel et al., 2003, ApJ, in press (astro-ph/0302209)) 
Silk, J., Rees, M.J., 1998, A&A, 331, LI (SR98) 

Tremaine, S., Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S.M., Filippenko, 
A.V., Green, R., Grillmair, C, Ho, L.C., Kormendy, J., Lauer, T., Magorrian, J., 
Pinkney, J., Richstone, D., 2002, ApJ 574, 740 

Yu, Q., Tremaine, S., 2002, MNRAS 335, 965 



This preprint was prepared with the AAS IM^jX macros v5.0. 



